Scaling of the linear response in simple ageing systems without disorder 
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The time-dependent scaling of the thermoremanent and zero-field-cooled susceptiblities in ferro- 
magnetic spin systems undergoing ageing after a quench to a temperature at or below criticality 
is studied. A recent debate on their interpretation is resolved by showing that for systems with 
a short-ranged equilibrium spin-spin correlator and above their roughening temperature, the field- 
cooled susceptibility XFc(i) — Xo ~ where xo is related to the equilibrium magnetization and 
the exponent A is related to the time-dependent scaling of the interface width between ordered do- 
mains. The same effect also dominates the scaling of the zero-field-cooled susceptibility xzFc(t,s), 
but does not enter into the thermoremanent susceptibility pTRM(t, s). However, there may be large 
finite-time corrections to the scaling of pTRM(t, s) which are explicitly derived and may be needed in 
order to extract reliable ageing exponents. Consistency with the predictions of local scale invariance 
is confirmed in the Glauber-Ising and spherical models. 

PACS numbers: 05.20.-y, 64.60.Ht, 75.40.Gb 



I. INTRODUCTION 



The comprehension of the physics of ageing phenomena is a topic of much current interest. While ageing was 
originally studied in glassy systems, there are many conceptual problems which are conveniently studied first in 
simple ferromagnetic systems. In this paper, we shall study a ferromagnet with a critical temperature > 0, initially 
prepared in a fully disordered (infinite-temperature) state which is quenched at time t — to a temperature T < Tc- 
We shall consider throughout a dynamics with a non-conserved order parameter. Physically, the ageing process 
proceeds via the growth of correlated domains of size L{t) ^ t^/^ and the slow motion of the domain boundaries 
drives in turn the slow temporal evolution of macroscopic obscrvables, see |^ |^ for recent reviews. It has turned 
out that ageing phenomena are more fully revealed through the study of two-time quantities such as the two-time 
correlator C{t, s) and the two-time linear response function R{t, s) defined by 
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where (j){t) is the time-dependent order parameter, h{s) is the magnetic field conjugate to (/), t is referred to 
observation time and s will be called the waiting time. Causality implies that R{t, s) = for t < s. 

Ageing systems may display dynamical scaling in the long-time limit Q, |^, |^ . Specifically, consider the two-time 
functions in the ageing regime t 3> tmicro, s 3> imicro and t = t — s ^ imicro, where tmicro is some microscopic time. 
Then one has the scaling behaviour 
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where the scaling functions fc,R{x) have the following asymptotic behaviour for a; 00 
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and A(7 and A/? are the autocorrelation 0,0 and autoresponse @ exponents, respectively. In general, the exponents 
\c,Ft, and z will take different values for T < Tc and for T = Tc- In particular, z = 2 for T < Tc and a non-conserved 
order parameter. 

The values of the exponents a and b are collected in tableland depend on the equilibrium spin-spin correlator Coq 
as follows I7I. If Coq(r) - e-l''l/« with a finite one says that the system is of class S, while if Ccq(r) ^ |r| ^ ^ 
the system is said to be of class L, where 77 is a standard equilibrium critical exponent. We point out that for systems 
of class S, the result a — 1/z follows from the well-accepted intuitive picture that ageing effects come from the slow 
motion of the domains walls which separate the well-ordered domains in systems undergoing coarsening 0, |3i ■ 
A different value for a would invalidate this physical picture. 

Recently, it has been proposed that ageing systems might possess a larger dynamically generated space-time sym- 
metry than mere dynamical scaling Indeed, for any given value of z, infinitesimal local scale transformations with 
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a space-time-dependent rescaling factor 1 + e(t, r) can be constructed. In particular, the following explicit expression 
for the response function is obtained from the condition that R{t, s) transforms covariantly under the action of local 
scale transformations 0, IT(tI | 



R{t,s)^rU- (t-s)-'-'^ (4) 



and where v'q is a normalization constant. This prediction has been confirmed in several models, notably the kinetic 
Ising model with Glauber dynamics 0, and several variants of the exactly solvable spherical model [Tol ITsIl , 
as is reviewed in 9, 14] . Furthermore, predictions exist for the spatio-temporal response. A particular simple result 
is found for the case 2; = 2 IS, [l3 
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LM is a non-universal constant) and has been confirmed in the 2D and 3-D Glauber-Ising model quenched to T < Tc 
[lj| (and also in the spherical model for d > 2). Since the derivation of (jSJ depends on the Galilei invariance of 
the coarsening process, its verification is strong evidence in favour of Galilei invariance as a dynamically generated 
space-time symmetry of phase-ordering. 

One would obviously like to be able to test statements of such a general nature as widely as possible, but since there 
are so few exactly solvable and non-trivial model of ageing, quantities such as R{t, s) must in general be obtained 
numerically. This paper addresses an important conceptual point on the interpretation of the scaling behaviour of 
several commonly used observables related to R{t, s) and we shall try to clear up a debate which has arisen in recent 
years. 

Numerical information on R{t,s) can be obtained using a by now standard method devised by Barrat Hid] . One 
perturbs the system by a random magnetic field hi with zero mean hi — 0. We shall use a binary field hi = ±h below 
in our Glauber-Ising model simulation, but a gaussian random field is also possible. Instead of measuring R{t, s) 
directly, two common procedures run as follows. Either one quenches the system and turns on the magnetic field 
after the waiting time s has elapsed and then works with the zero- field-cooled susceptibility Xzfc(^7 s). Alternatively, 
one may also keep the random field until the waiting time s when it is turned off and then has the thermorcmanent 
susceptibility pTRuit, s). These are related to R{t, s) as follows 

XzFc(i,s) = X{t,s)^ / duR{t,u) = s~y^{t/s) (6) 



PTBM{t,s) = p{t,s)= / duR{t,u) = s-^fM{t/s) (7) 

Here, we straightforwardly used the scaling forms Q without paying attention to the conditions of validity of these 
and in particular did not pay any attention as to whether t — s ^ imicro holds true or not. As we shall show, however, 
careful consideration of these conditions is crucial in order to obtain valid scaling forms for the integrated responses 
XzFc(i, s) and pTnuit, s). 

On a discrete lattice A C Z'', the integrated responses (|6I7|I are obtained by measuring the time-dependent magne- 
tization 



"^^(gWO.,) (8) 

where |A| is the number of sites of the lattice A. Depending on whether one works in the zero-field-cooled or the 
thermorcmanent protocol, one obtains xzfc = MzFc/h or p ttim = MTYiu/h, respectively. 

In recent years, Corberi, Lippiello and Zannetti (CLZ) \Vn. [T§ | have studied in great detail ageing in simple 
ferromagnets such as the Glauber-Ising model. Most notably, and based on MC data on xzfc in between one and 
four dimensions, they have argued that for the phase-ordering kinetics (i.e. T < Tc) in the Glauber-Ising model, 
the exponent a takes a value different from the generally accepted value a = 1/z ~ 1/2. Their studies consider 
the two-time autocorrelation C(t, s) and the dynamic susceptibilities xi^T s), p{t, s) of the 0(n) vector model. They 
separate the autocorrelator into a 'stationary' and an 'ageing ' part C{t, s) = C("*)(i - s) + C(^g<=)(t, s) and similarly 
x(t, s) = x'^^^H^ ~ ■5) + X^^'^'^HtjS), where x'^^'H^i •s) is defined such as to satisfy the fluctuation-dissipation theorem 
with C'''*^(t, s) (and similarly for p{t,s)). Furthermore, CLZ admit dynamical scaling for the 'ageing' part only and 
further introduce two distinct exponents and Op according to 

X^^^''Ht,s)^s-^-f^{t/s) , p(^^^\t,s)^s-^^Ut/s) (9) 
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If a„ is the value of a-^ in the 0(n)-model, CLZ propose 

_({d-l)/A ; d<3 _Ud^2)/2 ; d<4 

"i-\l/2 ; d>3 ' ; d>4 ^^^^ 

for the Glauber-Ising {n = 1) and spherical model {n = oo), respectively (additional logarithmic factors may occur 
at d = 3,n = 1 and (i = 4,n = oo, respectively). In order to account for these results, CLZ invoke a dangerous 
irrelevant variable and claim that d = 3 were a critical dimension of coarsening in the Glauber-Ising model. On the 
other hand, from a log-log plot of p(^sc) (j^^ g'j against s in the Glauber-Ising model ". . . no statement on Op can be 
made ..." 0, but roughly Op « 0.55 — 0.60 in 2D and Op « 0.65 — 0.70 in 3D, see 0, Fig. 6]. These conclusions 
are supplemented by a discussion of the necessary condition s ^ tmicro needed for the validity of the scaling Q and 
about the microscopic times <micro which may enter there „ISj. However, the crucial role of the additional condition 
t — s ^ tmicro which must also be satisfied for eq. (|2Jl to hold is not addressed by CLZ. 

Already from a pure phemomenological point of view and given the scaling relations |(5J and the definitions of x(t, s) 
and p(t, s), the assertion of CLZ of having distinct exponents ^ ap cannot be maintained. Rather, it is necessarry 
to go beyond a pure scaling analysis in order to understand how the ageing behaviour can be correctly extracted from 
integrated response functions. 

Our argument runs along the following lines. 

1. We consider the field-cooled susceptibility 

XFc(i) = XzFc(^, s) + PTRM(i, s) = / duR{t,u) (11) 

Jo 

For ferromagnets of class S, one has Xfc(0 ~ Xo ~ t~^j as we shall show in section Ull Here T xo — ^ ~ "^oq 
given by the mean equilibrium magnetization. The exponent A is a new exponent without a direct relationship 
to ageing, rather it is related to the roughness of the interface between ordered domains. We stress that since 
the power-law behaviour of xfc (t) is independent of the waiting time s it has no relation with a possible ageing 
behaviour. In particular, we find A — 1/4 for the 2D Glauber-Ising model with T < Tc- Furthermore, for a fixed 
scaling variable x = t/s one has, with the scaling function g{x) ^ x~'^ 

XZFc(i, s) = XFc(t) - PTRuit, s) - xo + s-^g{x) + 0(5"'') (12) 

Here we anticipate an important result of section Hvl namely that Ptrm(^is) ~ see also eq. (|13|l below. 
Since for rough interfaces one has A — a — A— l/z<0, it follows that ® cannot be used. Furthermore, 
the splitting of x{t,s) into a 'stationary' and an 'ageing' part advocated by CLZ is in contradiction with H12|l . 
Rather, the leading scaling of xfc with the waiting time s is unrelated to the ageing behaviour of the model 
as described by (0). Indeed, the terms describing ageing only occur as subleading terms in xzFc(i, s) and are 
therefore difficult to extract. In particular, the exponent A cannot be identified with the ageing exponent a. An 
example is provided by the 2D Glauber-Ising model, where a = l/z = 1/2 but A = 1/4. 

2. For ferromagnets of class L, we shall show that A = 0. We shall explicitly test this in section Hill for the non- 
equilibrium critical dynamics (i.e. T — Tc) in the 2D Glauber-Ising model and in the kinetic mean spherical 
model for any T < Tc. 

3. Eq. H12|l shows that ageing effects merely provide a finite-time correction to the scaling of xzfc- On the other 
hand, ageing terms are leading in the thermoremanent magnetization which scales more precisely as |^ 

PTRuit, s) = s-'^fMit/s) + s-^'^'^gMit/s) (13) 

where the scaling functions fnix) and guix) are related to the response function R{t,s) and can be found 
explicitly, see section IIVI For example, for a system in class S with an uncorrelated initial state, one has 
= ^ d/2 ^9J and since a = l/z, these two terms may be of almost the same order and a simple log-log 
plot may not be sufficient to yield a precise value of a for times accessible in present simulations. Indeed, this 
situation occurs in the 2D and 3D Glauber-Ising model. However, subtracting the leading correction according 
to eq. (|13|l allows to reliably determine a and the scaling function fuix) 0, We shall describe this in 

section IIVI 

In an appendix possible implications for the scaling of i?(t, t ~ e) for e — > are discussed. 
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II. SCALING OF XFC FOR SHORT-RANGED EQUILIBRIUM CORRELATORS 



We begin by discussing the time-dependent scahng for the field-cooled susceptibility XFc{t) for the Glauber-Ising 
model. We consider the Ising model on a hypercubic lattice with periodic boundary conditions and the equilibrium 
Hamiltonian Ti = — j) '^i'^j where the sum is over nearest neighbours only. We use heat-bath dynamics defined 
through the stochastic rule 

ai{t + 1) ±1 with probability i [1 ± tanh {Hi{t)/T)] (14) 

with the local field Hi{t) = X]y(i) '^y(0 where y(i) runs over the nearest neighbours of the sites i. Initially, the 
system is prepared in an infinite-temperature state. Thermoremanent and zero-field-cooled susceptibilities PTRuit, s) 
and XzFc(i7'S) can now be measured by perturbing the model by a binary random field hi = ±h with zero mean 
hi = and using eq. (jSJ. From these two independent measurements, we obtain the field-cooled susceptibility 

XFC (t) = PTRM {t, s) + XZFC {t, s) (15) 

In figure ^ we show xfc (t) so obtained for several values of the waiting time s in the case of the two-dimensional 
Ising model, at the fixed temperature T = 1.5 < Tc. Clearly, xfc is independent of s, which means that we are well 
inside the linear-response regime. Furthermore, within the times considered, our data are consistent with the scaling 
behaviour 

XFc{t)^t-^ , A « 0.25 (16) 

Since this dynamical scaling is completely independent of the waiting time s, it is unrelated to any ageing behaviour 
which might occur in this model. 

We now give a heuristic argument in order to understand where the scaling form and the value of A come from. 
Consider a simple ferromagnet in d > 1 dimensions which is quenched at time t — from an infinite-temperature 
initial state to a final temperature T < T^. The dynamics is assumed to be purely relaxational, i.e. without any 
conservation law. Microscopically, it is well-known that the configurations of the system consist of fully ordered 
domains of spins, of a typical size L{t) - with z = 2. We now perturb with a random field of zero mean hi — 
and wish to obtain the susceptibility xfc = Mpc/h from eq. First consider the case when T = 0. Then, because 
the spins deep inside the cluster are ordered, the only non-vanishing contribution to xfc comes from the spins from 
near the interfaces between the ordered clusters. We denote the interface density by pi{t) and have pi{t) ^ L{t)~^, 
see If w{t) is the interface width, we have 




XFcit) = r7n:^{l^<Jiit)hi) = j^{ ^ ai{t)hi) ^ L{t)-'w{t) (17) 

i(E interfaces / 

For a finite temperature T > 0, the order deep inside the clusters is not perfect and there remains a residual 
contribution to the susceptiblity. We then have, for large times 

XFc{t)^Xo + m-^w{t) , Txo = l-m,\ (18) 

where mcq is the equilibrium magnetization (for the 2D Ising model, nicq = (1 — sinh(2/r)^'*)^/^, see e.g. ^3)- The 
dynamics of a (d— l)-dimensional interface in a d-dimensional system (d > 2) can be described by the dynamics of a 
height model i21i] of continuous height variables S M and the equilibrium Hamiltonian Hlv] = — j') (''^j ~ ^j') 
with nearest-neighbour interactions and where r is the effective interfacial tension |21l. In adopting this description, 
we tacitly assume that the system is above its 'roughening temperature' Tr, see e.g. 1231, such that the fluctuations 
of the interface are unbounded (rough interface) for T > Tn and with bounded fluctuations for T < (smooth 
interface). This condition is always satisfied for d — 2, since then Tr = and indeed the description adopted here 
can be derived rigorously [l^. On the other hand, for d = 3 one has Tji w 0.5 Tc and finally, Tr — oo ior d > A. If 
the dynamics of the model is described by a Langevin equation, the squared interface width was shown by Abraham 
and Upton to scale for large times as w{ty = {vo{t)'^) ^ t^/'^ in 2D and w(i)^ ~ Ini in 3-D 21]. Inserting this into 
eq. (|17|l . wc find the following leading dynamical scaling for xfc in the Glauber-Ising model as t — > oo 

(i\ / ; if d = 2 

XFc(t)-Xo -^1,-1/2(1^,)!/. . if (19) 
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and provided that T > Tr. For dimensions d > A, and generically if T < Tr, one expects a flat interface with 
w{t) ^ const, and consequently XFcit) ~ Xo ^ t^^^^. Reconsidering figure^ it is easy to check that the final approch 
of xfc towards xo only occurs at times much beyond those accessed by our simulation. 

Combining eqs. (|12|l and H19|) . we have therefore reproduced the findings of CLZ in the Glaubcr-Ising model and 
have also made clear the physical origin of these results. Besides, a definite prediction for the logarithmic factor in 
d — 3 was obtained. In summary, we have for the Glauber-Ising model A = 1/4 in 2D and A = 1/2 for all d > 3, up 
to a known logarithmic correction in 3-D. 



III. SCALING OF XFC FOR LONG-RANGED EQUILIBRIUM CORRELATORS 

We now ask whether the heuristic discussion of the scaling of xfc presented in the previous section can be taken 
over for systems of class L. Indeed, the main physical difference with respect to systems of class S is that although 
correlated clusters of size L{t) ~ t^^^ form, fluctuations do persist in the interior of these clusters on all length scales 
up to L{t). This means that one should consider an 'interface width' scaling as w(t) L{t). This in turn leads to 
Xfc(^) const, and A = Q (on the other hand, since the clusters should have no 'inside', we do not expect a term xo 
to occur). 

We now test this heuristic idea in the exactly solvable mean spherical model and shall also present evidence from 
the 2D Glaubcr-Ising model quenched onto criticality. 

1. First, we consider the mean kinetic spherical model, see e.g. 0, 0) 13 113 ■ To each site x of a hypercubic 
lattice K dU^ one attaches a continuous spin variable S'x(i) S K, subject to the mean spherical constraint 

E (^x) = 1 (20) 
xeA 

By analogy with the Glauber-Ising model, we also add a random magnetic field /ix and the equations of motion read 

^5x(i) - E W - (2d - l{t))S^{t) + h^ + r,^(i) (21) 
y(x) 

where the noise and the random field have zero average (yyx) = = and the correlators 



{v^{t)Tiy{t')) ^2T6^^y5{t-t') , /ix/iy = 2r5x,y (22) 

where the temperature T and the width F are constants and i{t) is a Lagrange multiplier to be determined below. 
In addition, the noise and the field are assumed independent, i.e. ('7x(0'*y) — 0; ^^^^ i^^ addition the initial state is 
uncorrelated in the sense that (S'x(0)/iy) = 0. Here and in the following the average (X) is always taken over the 
initial conditions and the noise, while the average X is over the random field. The solution of this model follows 
standard lines [l^. Taking Fourier transforms, the solution of eq. (|21|l is 



5(q,t) 



exp(-w(q)i) 



5(q, 0) + fdt' e-('i)*' (/.(q) + ^^(q, t') 







(23) 



with the dispersion relation Ci;(q) — 2^^^-^(l — cosq^) and g{t) ~ exp(2 At' lit'))- The spin-spin correlator 



C(q, q'; s) = t)S{q' , s)j is readily found and we have in direct space the autocorrelator 

C{t, s) = Cx.x(t, s) = -i^ / dqdq' e'i^+'i'y- C(q, q'; t, s) 
(271-)-^" Jb2 

^ - gin) + 2T j\t' j\\s' / {^-±1^1-^^ 



A I tLl ]+2T I duf 



1 

~ V9{t)9is) 

where B is the Brillouin zone and with the definitions 



(24) 



A{t) = (2^)--^ / dqe-2'"('i)*C(q,0) (25) 
fit) = {27r)-^ I dqe-^'^f'i)* = {e^'* 10^)" (26) 
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and where Iq is a modified Bessel function. For infinite-temperature initial conditions A(t) = f{t). The mean spherical 
constraint (|2U|) gives C{t,t) — 1 and this leads to the following generalized Volterra integral equation 



git) = A{t) + 2T fdu fit - u)giu) + 2T f du' f du" f(t~ ] ^giu')giu") 



(27) 



which determines git). Finally, the response function is given by the usual equation s) = fi{t— s)/2)\/qis)/ qjt). 
At zero temperature T = 0, eqs. (|24I27|I are identical to those found for the spherical spin-glass p3.l2g|. 
The field-cooled susceptibility is given by 



XFc(i) 



/ duRit,u) = [ duf 
Jq Jo 



t-u \ giu) 

2 J^git) 



(28) 



At this point, it is instructive to re-derive the equivalence between the definition H28|l and eq. originally proposed 
[T^ for the Glauber-Ising model, in the context of the mean kinetic spherical model. Indeed, we have (|A| denotes 
the number of sites of the lattice A) 



|A| 



\xeA 



-i— ^ / dqdq' e'('i+'i (^(q, <)Mq'; 

I dqdq'e'('i+'i')- f dt' e^^^^^^' -'hl^ -hi^M^) 
1^1 ^Jb^ Jo V 9^^> 



= (27r)-''2r / dt' I dqe' 
Jo Jb 

= 2rxFc(0 



V 9it) 



(29) 



as asserted, and where we used in the second line eq. (|23ll and in the third line the field correlator 1221) . 

After these preparations we can test our heuristic picture. Using the techniques described in t25j . we obtain git) 
by solving eq. H27|l numerically. In figure |21 we show xfc(0 for the three-dimensional case, starting from an infinite- 
temperature initial state. We clearly see that XFc(i) saturates rapidly. Consequently, ^ = for all temperatures 
T <Tc. Similar tests can be performed for other values of d as well. 

We pause a moment in order to discuss the functional form of the response function. If F = 0, solution of the 
spherical constraint gives for T <Tn the well-known exact result, valid for all values of d and in the ageing regime 
s> 1 andt-s> 1 is, see H E iH US 111 IH 



i?(i,s) 



. d/4 



-d/2 



(30) 



with t'q = (47r)~''/^. From (O, we read off a = id — 2)/2 and A^/z — d/A. We point out that this exact result 
for a is in contradiction with the claim pO|l raised by CLZ. Furthermore, we see that the exact result eq. (E2J) has 
precisely the form Q predicted by local scale invariance 0| . A similar test can be performed for T = 0, ^3 ' 
or even in a spherical model with spatially long-ranged interactions It is important to note that the local scale 
invariance prediction applies to the full response function i?(t, s) and not to the part remaining after subtraction 
of a 'stationary' term, as suggested in p^. 

2. As a second example, we give numerical evidence that A = in the 2D Glauber-Ising model quenched to its 
critical point T = Tc- We used a standard heat-bath algorithm and measured the integrated linear response through 
eq. (jHJ. Data for xfc(^) thus obtained are displayed in figure |2I and we see that saturation occurs. 

In summary, extending a heuristic argument of sectional we have argued that for ferromagnetic systems in class L, 
one should find saturation for the field-cooled susceptibility, viz. xfc(0 ^ 0(1). We have confirmed this expectation 
in some models. 

The heuristic arguments in this and the preceeding section can only be applied in d > 2 dimensions. Indeed, the 
ageing behaviour of the ID Glauber-Ising model at its critical point Tc = is peculiar and will be discussed in the 
appendix. 
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IV. SCALING OF THE THERMOREMANENT RESPONSE 



Since we have seen that the terms which describe the physicahy interesting ageing effects are only sub-leading in the 
zero-field-cooled susceptibility Xzfc(^jS), we discuss in this section the scaling of the thermoremanent susceptibility 
PTR iuf^i s). As we shall see, the case when the waiting time s is small needs particular consideration. It is well-known 
|29| that the response with respect to a fluctuation in the initial state scales as 

R{t,0) ^t-^'^/'' (31) 

On the other hand, it can be shown that there is a time-scale tp ~ s'' with < C < 1 such that if the time 
difference t ~ t — s < tp, then the response function is still the one of the equilibrium system R(t, s) ~ i?gg(r) while 
scaling sets in and lO holds if t > tp. For example, in the spherical model one has C = 4/((i-|- 2) |23|. We then have 

PTRM{t,s) = / duR{t,u)^ / dr R{t,s ~ t) 
Jo Jo 

= / dTRgg{2T)+s dvR{t,s{l-v))+ dTR{t,S-T) 

Jq Jt„/s Jt, 



tp/ 

dr i?gg(2T) + ^ di; (1 - v)-^-Vr (^^j + Coo t"^«/^ (32) 

Here, following 0, ESli have introduced a third time-scale t,, such that s — = 0(1). In the third line, the limit 
s ^ oo is taken and the last term is estimated from the mean- value theorem, where Cco is a constant. Since we are 
perturbing with a random magnetic field, the system is not driven to a new equilibrium state and consequently, the 
first term poo = /g dTi?gg(2T) vanishes if the initial mean magnetization was taken to be zero. As a result we have 
the scaling form, already announced in 

PTRuit, s) = s-''fM{t/s) + s-^'^'^-gMit/s) (33) 

where fnix) and gM{x) are scaling functions. These results only depend on the assumption of dynamical scaling. If 
in addition local scale invariance applies, the form of the response function i?(i, s) is given by eq. and therefore 

Hi 







\ Z Z X J 



gM{x) ~ rix-^«/^ (34) 



and where ro,i are non-universal constants and 2F1 is a hypergeometric function p8|. 

The importance of taking this finite-time correction into account is illustrated in figure 0] For a fixed value of 
X = t/ s and T < Tc, we plot data for the Glauber-Ising model in 2D and 3D, respectively and the exact solution of 
the Langevin equation of the spherical model, with F = 0. A fully disordered initial state was used. We compare the 
data with the leading scaling form p{t, s) ~ s~° which for the times accessible does not fully describe the data, but 
inclusion of the second term in (|33|l gives a very good fit. From this fit we find the non-universal values of ro.i listed in 
tabled In order to achieve this, however, a must take the values given in table Specifically, for the Glauber-Ising 
model with d> 2 (which is in class S), we must take a = 1/2, while if we had chosen a = 1/4 as advocated by CLZ, 
only a fit of very low quality is obtained, see 0, H^l ■ This provides further evidence against the proposed eq. ifTUIl . 
We point out that having fixed ro,i for a given value of x, the scaling of ptrm(^, s) is completely fixed and does not 
contain any further tunable parameter. Repeating the comparison between the data and (|33|l for other values of x 
therefore allows to test the predictions of local scale invariance and this has been carried out in detail for the 2D 
and 3D Glauber-Ising model 0, 0| • While in the 2D Glauber-Ising model the curvature of the data, see figure 
already suggests the presence of strong finite-time corrections to scaling, in 3D it might appear at first sight that a 
reliable value of the exponent a could be derived, see figure^. However, the inclusion of the leading correction term 
from eq. H33|) together with the study of the x-dependence of the coefficients, as carried out in 7], is required in order 
to avoid a systematic error in the determination of a. In this way the expected value a = l/2 = l/2is reconfirmed. 
Carrying out the analysis of the scaling functions as described in j2i] , we can conclude that local scale invariance holds 
in the 2D and 3D Glauber-Ising model with T <Tc- 

A closer inspection of tableElshows that in the spherical model with d > 4 and T < Tc, the signs of ro and ri change 
and that the term which was treated as a correction in (|33|l becomes dominant. Indeed, the exactly known exponents 
a = {d — 2)/2 and Xr/z = d/A cross at the upper critical dimension d* = 4. We conclude that the scaling argument 
used in to derive the relation a = {d — 2 + ri)/z for class L systems should rather be viewed as a hyperscaling 
argument. 
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V. CONCLUSIONS 

In this paper, we have analysed the long-time behaviour of integrated response functions in ageing spin systems 
without disorder. While from the scaling form (|2Jl, one might have expected a simple scaling x(i, s) ~ s"" of the 
integrated response, we have shown that matters are more complicated. Conceptually, the issue can be clarified by 
studying the scaling of the field-cooled susceptibility XFc{t) and we have seen that two broad classes of systems must 
be distinguished, which are called classes S and L, according to whether their equilibrium spin-spin correlator shows 
short-ranged or long-ranged spatial decay, respectively ^7J . Specifically, we have found 

1. For systems of class S, we have 

XFc{t) - xo - t-^ (35) 

where Txa = 1 — ''^cq and rrieq is the equilibrium magnetization.The exponent A is related to the interface width 
exponent k. If u'(t) ~ i", we have seen that 

A = ^-K ; (class S) (36) 

but provided that the temperature T > is above the roughening temperature T/j. 

2. For systems of class L, eq. still holds with xo = 0, but with 

^ = ; (class L) (37) 

3. The scaling of the thermoremanent susceptibility PTRM(i, s) may be affected by a relatively large finite-time 
correction term, see eq. (|33f) where the associated scaling functions can be found explicitly from local scale 
invar iance. 

4. Consequently, the leading time-dependence of the field-cooled susceptiblity xzFc{t, s) ~ Xo ~ Xfc(^) ~ Xo ~ 
whereas the ageing terms expected from integrating eq. l(2Jl merely arise as finite-time corrections. 

Our results reproduce the entirety of known results in ageing ferromagnetic spin systems with d > 2 but do not 
require to postulate an upper critical dimension d = 3 in the Clauber-Ising model. For systems of class S, the 
exponent relation a = 1/z is reconfirmed while the proposed eq. (|1U|) is invalidated. 

If questions of simulational efficiency play no role, it might be technically easier to avoid both the ZFC and the 
TRM protocol, as already suggested in [l^. For example, we propose the 'intermediate' protocol, which runs as 
follows: quench the system at t — without a magnetic field and fix a waiting time s. At time s/2, turn on a random 
magnetic field and keep in on until the waiting time s. Then turn the field off again and measure the magnetization 
at the observation time t > s. The intermediate integrated response is 

piNT(t, s) = AfiNT(t, s)/h du R{t, u) = s-'^/iNT(i/s) (l + o(s-^«/^)) (38) 

and should be free of the leading term coming from the interface roughness as well as the finite-time correction of order 
0{s~^"-/ '^). We illustrate this for the mean spherical model in figureEl where it can be seen that already for times 
much shorter than in figure^, the linear response ()28|l obtained from the exact solution of the Langevin equation 
(with F = 0) converges to the expected power law, with a = 0.5 in 31?. It would be interesting to see if a recent 
method to calculate xzfc(^j s) directly in the Glauber-Ising model 31] (which in turn is based on a method to estimate 
R{t, s) in an Ising model with a different dynamics modified from Glauber dynamics H^) could be generalized to find 
PiNT(i, s) as well. 
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Appendix 

In the main text, we studied the time-dependent scaling of the integrated response XFc{t). We shall now attempt 
to extract further information on the scaling of _R(t, s) when t k, s. Note that the scaling form eq. Q cannot be used 
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close to the upper integration limit, since there the condition t = t — s ^ 1 needed for the validity of |5J does not 
hold. Rather, one might argue that there should be some time-scale t* such that j^Jl holds for waiting times s <t*. 
If in addition t — t* is some 'small' constant At* for large times, we may write 



XFc{t) = ^ dsR{t,s) + J^dsR{t,s) 



dvv''-^fR{v) + {t-t*)R{t,t') 

Jt/t* 

poo 

r"" dw""-^ fBiv) + At* ■ R{t,t') (Al) 



where we used the mean- value theorem and t' E [t* , t]. If the second term in I|A1|I is the leading one, then comparison 
with the scaling of xfc gives R(t,t) ~ 1 -I- 0{t~"^), provided the conditions stated above are satisfied. Here the 
value of the constant is fixed from the physical consideration that for t — t — s <^ 1 the system should still be in 
quasiequilibrium. 

We can now try and see to what extent this argument applies in specific models and to what extent it might be 
justified. The only model of class S we studied here is the 2D Glauber-Ising model at T < Tc. Then a — 1/2 and 
A = 1/4 and the contributions close to the upper integration limit are indeed dominant. One should therefore expect 
lime^o R{t, t-e)~ cste. + 0{t~^/^) for t large. 

On the other hand, we have considered in section llTTl several examples of systems of class L, where w e have found A = 

throughout. First, for the spherical model, the exact expression for R{t, s) gives R{t, 2e) = f{s)\/g{t — 2e) I g{t) ~ 

1 -I- e (/'(O) — g'{t)/g(t)) + O(e^). Since for t large, one has g(t) ~ t"^ with a known exponent ■0 we obtain for 
e — > R{t, t — 2£) ~ cste.-|-0(i~^), and have indeed A — Q, since the constant is different from unity. Second, a similar 
conclusion should also hold true for the 2D Glauber-Ising model at criticality. Third, an analogous result can also be 
obtained in the 2D linear voter model at criticality ^33j . where the exact response function lim^^o R{t+£, t) {lnt)~^. 

Additional studies to test these conclusions further would be welcome. 

The heuristic arguments presented in this paper assume d > 2 space dimensions throughout. As an illustration 
what can happen in one dimension, we now consider the ID Glauber-Ising model at its critical temperature = 0. 
This model has been investigated countless times and we shall not repeat the details of its definition (see e.g. 
and references therein) here. The exact autoresponse function is 

R{t, s) = e-(*-^)/o(t - s) (1 - Ci(s)) (A2) 

where Ci (s) is the equal-time correlator of two spins on neighbouring sites. For an infinite-temperature initial state, 
one has ISa 



CAs)=e 



-2s 



h{2s) + 2j2ln+i{2s 

n=l 



e-^q/o(2s) + /i(2s)] (A3) 



where in the second step the Bessel-function identity — Io{z) + 2J2'^=i In{z) was used. It has recently been shown 
that for any initial correlators of the power-law form Cini(f') '^-^ r"'^ with z/ > the leading scaling behaviour of Ci(s) 
for a large is unchanged [35l |. Therefore, for any power-law correlated initial state we have for t large 

\\mR{t,t~e)~-^-t-^'^ (A4) 

Comparison of eqs. I|A1(I and l|A4|l . if admissible, would then give A = 1/2. Such a result would appear to be natural 
in the setting of section |nl if we recall that at T = there are fully ordered domains where the width of the domain 
walls is one lattice constant, thus w{t) — cste. 

On the other hand, integrating (|A2IA3|I . the field-cooled susceptibility becomes for an uncorrelated initial state 

XFc(0= / dse-'-'h{t-s){h{2s)-Vh{2s)) (A5) 
Jo 

and we show this in figure Clearly, there is saturation for large times and in the spirit of section IIIII one would 
conclude A — Q. Therefore, a naive application of the arguments valid for d>2 would lead to two different values of 
A in this ID model. We leave open the question how this apparent inconsistency might be resolved. 
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TABLE I: Non-equilibrium exponents a and 6 for quenches from a fully disordered state onto and below the critical point 
Tc > of simple ferromagnets of the classes S and L defined in the text, according to 0- 







a 
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Class 


T = Tc 
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n)/z {d- 


2 + 7))/z 
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T <Tc 
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TABLE II: Values of the autoresponse exponent Ah and of the parameters ro and ri in the Glauber-Ising model in two 
dimensions at T = 1.5 and in three dimensions at T = 3 and of the mean spherical model at T = 2 in several dimensions. An 
infinite-temperature initial state was used. The dynamic exponent z — 2. 



model 


d 




ro 


ri 


Glauber-Ising 
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1.26 


1.76 ±0.03 


-1.84 ±0.03 
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1.60 


0.10 ±0.01 


0.20 ±0.01 


spherical 
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1.50 


0.180 ±0.01 


-0.081 ± 0.002 




3.5 


1.75 


0.20 ±0.01 


-0.126 ±0.003 




4.5 


2.25 


-0.056 ± 0.003 


0.095 ± 0.002 
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FIG. 1: Field-cooled susceptibility XFc(i) ~ PTRmit, s)+xzFc(i, s) in the 2D Glauber-Ising model at temperature T — 1.5 < Tc 
and for several waiting times. An alternating field hi — ±/i with h = 0.05 was used. The power law XFc{t) ~ t''^'^^ is also 
shown for comparison. The simulated systems contained 300 x 300 spins, and the data have been obtained after averaging over 
at least 5000 different runs. 




FIG. 2: 

Field-cooled susceptibility xfc (t) for the 3-D kinetic mean spherical model in a gaussian random magnetic field of width F 
and at r = 2 < Tc (full curve) and at T = ~ 3.96 (dashed curve). 



= 0.01 
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FIG. 3: Field-cooled susceptibility XPc{t) in the 2D Glauber-Ising model at the critical temperature T = Tc ~ 2.269. An 
alternating field hi = ±h with h = 0.05 was used. Systems with 300 x 300 spins have been simulated. 




FIG. 4: Scaling of the thermoremanent susceptibility p(t, s) = MTRM{t, s)/h as a function of the waiting time s for a fixed value 
of a: = t/s for three models, (a) The 2D Glauber-Ising model at T = 1.5 and with x = 7 (points) is compared with the scaling 
prediction eq. 1331 (full grey curve) and with the simple scaling form p{t,s) ~ s~^^^ (dashed curve), (b) Same comparison for 
the 3D Glauber-Ising model at T = 3 and with a; = 5. (c) The "AD spherical model at T = 2 and x = 5. The full curve is 
the exact solution of the Langevin equation, the dash-dotted line the scaling prediction eq. (1331 and the dashed line the simple 
form p{t,s) ^ s~^^^. 
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FIG. 6: Field-cooled susceptibility xfc (i) for the ID Glauber-Ising model at T = 0. 



